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Abstract 

We study the optimal transmission scheme that maximizes the local capacity in two- 
dimensional (2D) wireless networks. Local capacity is defined as the average information 
rate received by a node randomly located in the network. Using analysis based on analytical 
and numerical methods, we show that maximum local capacity can be obtained if simultaneous 
emitters are positioned in a grid pattern based on equilateral triangles. We also compare this 
maximum local capacity with the local capacity of slotted ALOHA scheme and our results 
show that slotted ALOHA can achieve at least half of the maximum local capacity in wireless 
*0 ' networks. 
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1 Introduction 

Seminal work of Gupta & Kumar pQ and the following works, e.g., [21 [3] quantify the capacity 
in wireless networks in the form of asymptotic scaling laws. However, these results may not be 
very useful for network protocol designers in comparing different medium access schemes that have 
different protocol overhead but follow the same scaling behavior. Our goal is to investigate the 
medium access scheme which optimizes the local capacity. Note that any such scheme may have 
no practical implementation but its evaluation is interesting in order to establish an upper bound 
on the local capacity in wireless networks. In our analysis, we will use first and second order 
differentiation of local capacity to prove that simultaneous emitters arranged in a grid pattern are 
locally optimal and, in 2D wireless networks, only square, hexagonal and triangular grid patterns 
are most optimal patterns. 

This article is organized as follows. In Section |2l we will discuss the model of our wireless network 
and define the local capacity. We will summarize some important related works in Section [21 The 
optimality of grid pattern based medium access schemes will be discussed in Section H] and their 
local capacity will be analyzed in Section [51 Section [6] will discuss the local capacity of simple 
ALOHA based scheme. In Section [71 we will evaluate the local capacity of grid pattern schemes 
and slotted ALOHA and concluding remarks can be found in Section [HI 



*INRIA Paris-Rocquencourt, France. Email: salmcin.malik@inria.fr 
tlNRIA Paris-Rocquencourt, France. Email: philippe.jacquet@inria.fr 



1 



2 System Model 



We consider a wireless network where nodes are uniformly distributed over an infinite plane centered 
at origin (0, 0). We assume that time is slotted and at any given slot, simultaneous emitters in the 
network are distributed like a set of points, S = {z\, z 2 , . . . , z n , . . .}, where Zi is the location of 
emitter i. The distribution of set S depends on the medium access scheme employed by the nodes 
and we only assume that, in all slots, the set S has a homogeneous density equal to A. 

Let Pi be the transmit power of node % and 7^ be the channel gain from node i to node j such that 
the received power at node j is Pciij. Therefore, transmission from node % to node j is successful if 
the following condition is satisfied 



^0 + E/^ Pkikj 

where N is the background noise power and (3 is the minimum signal to interference ratio (SIR) 
required for successfully receiving the packet. We assume that all nodes use unit nominal transmit 
power and we only consider large-scale pathloss characteristics, i.e., 7^ = \zi — Zj\~ a , where a > 2 
is the pathloss exponent and |.| is the Euclidean norm of the vector. We also assume that Nq is 
negligible. Therefore, the SIR of emitter i at any point z in the plane is given by 



Si(z) = 





z - 


- Zi\ 


—a 




z — 


Zj\ 


—a 



(1) 



We call the reception area of emitter i, the area of the plane, Ai(\,(3,a), where this emitter is 
received with SIR at least equal to (3. The area, Ai(X, j3,a), also contains the point Z{ since here 
the SIR is infinite. The average size of Aj(A, 0, a) is a(X, (3, a) and it is independent of the location 
of Zi. 

We are interested in local capacity which is defined as the average information rate received by a 
node randomly located in the network. Consider a node at a random location z in the plane and 
let N(z, (3, a) be the number of reception areas it belongs to. The expected value of N(z, (3, a) is 
given by [1] 

E(N(z,P,a)) = \a(\,/3,a) . (2) 

The average information rate received by the node, c(z,(3,a), is equal to E(N(z, j3, a)) multiplied 
by the nominal capacity. We assume unit nominal capacity and we have 

c(z, f3, a) = E(N(z, P, a)) = Aa(A, p, a) . (3) 



3 Motivation and Related Works 

In related works, focus has been on the medium access schemes like ALOHA, carrier sense multiple 
access (CSMA) or, in some instances, node coloring as well. Some of these works are as follows. [5] 
studied slotted ALOHA using a very simple geometric propagation model. Under a similar propa- 
gation model, [6] evaluated CSMA and compared it with slotted ALOHA in terms of throughput. [7] 



2 



used simulations to analyze CSMA under a realistic SIR based interference model and compared it 
with ALOHA (both slotted and un-slotted). [8HT0] introduced the concept of transmission capacity, 
defined as the maximum number of successful transmissions per unit area at a specified outage prob- 
ability, and evaluated ALOHA, CSMA and code division multiple access (CDMA) protocols. [TT] 
analyzed local (single-hop) throughput and capacity with slotted ALOHA, in networks with random 
and deterministic node placement, and TDMA, in ID line-networks only. [12] determined the opti- 
mum transmission range under the assumption that interferers are distributed according to Poisson 
point process whereas [13J gave a detailed analysis on the optimal probability of transmission for 
ALOHA which optimizes the product of simultaneously successful transmissions per unit of space 
by the average range of each transmission. In contrast to these works, we investigate the most op- 
timal medium access scheme which optimizes the local capacity in wireless networks. We will also 
compare this optimal scheme with slotted ALOHA. More detailed comparison with other schemes 
will be done in the continuation of our work. 



4 Optimality of Grid Pattern Schemes 

It can be argued that optimal local capacity in wireless networks can be achieved if simultaneous 
emitters are positioned in a grid pattern. However, designing a medium access scheme, which ensures 
that simultaneous emitters are positioned in a grid pattern, is very difficult because of the limitations 
introduced by wave propagation characteristics and actual node distribution. Specification of a 
distributed protocol that allows grid pattern transmissions is beyond the scope of this article. Note 
that wireless networks of grid topologies are studied in, e.g., [HJ[T5]. I n contrast, we assume that 
only the simultaneous emitters form a regular grid pattern. 



4.1 Definitions 

In order to simplify our analysis, we define the following functions. 

1. The density of the set S is given by the limit as 

u(S) = lim h z \<R ■ 

V ' R->oo 7TR 2 ^ |2sl - 
i 

2. We define a function g(z) as 

9i(z) 



where a > 2. The function gi(z) is similar to the SIR function Si(z) except that the summa- 
tion in the denominator factor also includes the numerator factor. In order to simplify the 
notations, we will remove the reference to z when no ambiguity is possible. 
We define a function /(gi) as 

f{9i) = l 9i{z)>P' > 

for some given /3'. We consider without loss of generality that the value of /?' is given by 
j3' = -p-^. Therefore, if transmitter i is received successfully at location z (i.e. with SIR at 
least equal to 0), then gi(z) > (3' and /(<?«) is equal to 1. 
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4. Reception area of an arbitrary emitter i is given by 



ai = j f(9i)dz 2 . 

Note that the integration is over the plane and the notations are simplified by taking dxdy 
equal to dz 2 . 

5. We define a function h{z) as 

M*) = £/G&) , 

ies 

and it is equal to the number of emitters which can be successfully received at z. Note that 
if ft > 1, the maximum value of h(z) is 1. 

6. We also define E(h(z)) by the limit as 



E(h(z)) = hm -L / h{z)dz 2 , 



with 



lim — i 



E(cr,) = lim — > a 



where n is the number of emitters on a disk of radius R centered at (0, 0). As R approaches 
infinity (because of an infinite plane), n approaches infinity. Our objective is to find the 
spatial distribution of the set S which optimizes 'E(h(z)). Note that E(/i(z)) is equivalent to 
E(N(z, ft, a)) and local capacity given by expressions ([2]) and ([3]) respectively. 

4.2 First Order Differentiation 

Let us denote the operator of differentiation w.r.t. Z{ by V«. For i ^ j, we have 

Z - Zi 



z — z, 



|2 



and 



Therefore 



Vtfi = a(g 2 -gi)y- — - 

\Z - Zi 



2 



Vih(z) = V, f(9i) = f'(gi)Vi9i + f'(9j)Vi9j 
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We know that J h(z)dz 2 = oo, we nevertheless have a finite Vi / h(z)dz 2 . In other words, the sum 
^j^i converges for all z. For all j in S, Yli^i a j = 0- Indeed this would be the differentiation 
of <jj when all points in S are translated by the same vector. Similarly, 

h{z)dz 2 = . 

i J 

Theorem 4.1. If the points in the set S are arranged in a grid pattern then Vi J h(z)dz 2 = 
Vi a j = 0; an d grids patterns are locally optimal. 

Proof If S is a set of points arranged in a grid pattern, then: V« / h(z)dz 2 = ^ • VjCTj would be 
identical for all i and, therefore, would be null since ^ Vi f h{z)dz 2 = 0. 

We could erroneously conclude that, 

- all grid sets are optimal and 

- all grid sets give the same 'E(h(z)). 

In fact this is wrong: we could also conclude that E(o"i) does not vary but this will contradict that 
v{S) must vary. The reason of this error is that a grid set cannot be modified into another grid 
set with a uniformly bounded transformation, unless the two grid sets are translated by a simple 
vector. □ 



4.3 Numerical Differentiation of First and Second Order 

In order to prove that grid patterns are also locally maximum, we must show that V^ J h(z)dz 2 < 0. 
Obviously, analytical formulation to prove this property is very challenging. However, we can 
develop a numerical differentiation technique to show that this is true in case of grid patterns. 

We know that J h(z)dz 2 = oo, and we are only interested in the behavior of Vi J h(z)dz 2 and 
V^ J h(z)dz 2 . Therefore, we define a function U as 

U = [ h{z)dz 2 , (4) 

J A 

where the integration is over a large area, A, with location of emitter i, Zi, at the center of this 
area. Note that the area, A, is large enough so that a slight perturbation in the location of Zi may 
have infinitesimal effect on the reception areas of the emitters near the edges of A. Let Uij denote 
the result of the integration in (j3J) and we denote the differentiation of U w.r.t. Zi as 

v ' u -(^ u 'W. u )^ (u " u ' h 

We will compute VjC/ and VfU with slight perturbation in Zi and use the following method. 

Using central difference equations, we can write U x ~ 2Aa;(^+i,j — ^»-ij)> anc ^ U y ~ 2^(^m+i — 
Uij-i), where Ax and Ay are slight perturbations in the location of along x— axis and y— axis 
respectively. The value of C/j+ij is computed from (@| with perturbed position of emitter i given 
by z\ = (xi + Ax,yi), and this also applies to other notations. 
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Similarly, the second order partial derivatives can be written as U xx m -J-^(Ui + i t j — 2Uij + Ui-ij), 
U vv ~ -^(Uij+i-^Uij+Uij^), and mixed derivative can be written as U xy = U yx « A ^ xAy {U i+ld+x - 

In order to prove that grid patterns are locally maximum, we will analyze the eigenvalues of the 
2x2 Hessian Matrix, 



H{U) 



U X X U X y 
Uyx Uyy 



In other words, we can prove that the grid patterns are locally maximum if, in this case, we can 
show that the determinant, \H(U)\, is greater than zero and U xx < 0. 

We will introduce the analytical method to compute U and its first order and second order partial 
derivatives in £J5] and present the results of numerical differentiation in §7.11 

We prove that the grid sets are locally optimal and maximum within sets that can be uniformly 
transformed between each other. In order to cope with uniform transformation and to be able to 
transform a grid set to another grid set, we will introduce the linear group transformation. 



4.4 Linear Group Transformation 



Here, we assume that the points in the plane are modified according to a continuous linear transform 
M(t) where M(t) is a matrix with M(0) = I, e.g., M(t) = I + tA where A is a matrix. Without 
loss of generality, we only consider cq, i.e., the reception area of the emitter at z which can be 
located anywhere on the plane. Under these assumptions, we have 

^°"o = ^(A^.Vjcro) = tr(^ A T Zi <g> ViOo) . 

i i 

In other words, using the identity B ) = B, the derivative of o"o w.r.t. matrix A is exactly 



equal to D = £V Zi <g) VjCJo, such that D 
We can write the following identity 



Dxx D 

Dyx D 



yy 







tr(A^-, 



d_ 

dt 



<ro(t, A) 



t=0 



where <Jo(t, A) is the transformation of o"o under M(t), i.e., <Jo(t, A) = det(I + At)cr . We assume 
that M(t) = (1 + t)l with A = I, i.e., the linear transform is homothetic. 

Theorem 4.2. D is symmetric and tr(D) = 2<To. 

Proof. Under the given transform, a (t, A) = cr (t,l) = (1 + t) 2 a . As a first property, we have 
tr(D) = 2<7o, since the derivative of cr w.r.t. identity matrix I is exactly 2cq (i.e., tr(A T D) = 
tr(D) = <Xq(0,I) = 2o~q). The second property that D is a symmetric matrix is not obvious. The 

easiest proof of this property is to consider the derivative of ctq w.r.t. matrix J = 



1 



which 



is zero since J is the initial derivative for a rotation and reception area is invariant by rotation. 
Therefore, tr(J T D) = D yx — D xy = 0, which implies that D is symmetric. □ 



6 



Note that D can also be written in the following form 



D = E 



Vi<7o 



Zi ® Vi/(ff ) 



Let T be defined as T = / dz 2 - z { ) <g> Vif(g ) , such that D = / Ei ^ ® Vif(g )dz 2 - T . 
The purpose of these definitions will become evident from theorems 3 and 4. 

Theorem 4.3. We will show that f z <g> Vif(go)dz 2 is equal to a I and, therefore, D = o" 7— T. 
VKe will also prove that T is symmetric. 

Proof. From the definition of T, we can see that the sum J2i( z ~ z i) ® ^ifi.9o) leads to a symmetric 
matrix since 

,2 



a 



J \\z-Zq\ 2 



'z - Zq) ®(z- Zq) + 



E 



9o9i 



Z - Zi 



:(z - z^ ®(z- Zi)^dz 2 , 





O"0 










O"0 



and the left hand side is made of (z — Zi) <S> {z — Zj) which are symmetric matrices. This implies 
that T is also symmetric. 

We can see that £V ^if(9o) = ~^f(.9o), and using integration by parts we have 
J ^z®Vif{g Q )dz 2 = -lx 

i 

J x£f(g )dxdy J x-§-f(g )dxdy 
Jy£f(9o)dxdy J y£f(g )dxdy 

which is symmetric and equal to er I. The sum/difference of symmetric matrices is also a symmetric 
matrix and, therefore, D is a symmetric matrix and D = <7oI — T. □ 

Now, we will only consider grid patterns and, by virtue of a grid pattern, we can have 

E(<7i) = a = J f(g )dz 2 , 

and E(h(z)) = v{S)o~o. Under homothetic transformation, v{S) and o"o are transformed but v{S)(Jq 
remains invariant. 

Theorem 4.4. If the pattern of the points in set S is optimal w.r.t. linear transformation of the 
set, D = a I and T = 0. 

Proof. The derivative of a w.r.t. matrix A is exactly equal to D. Similarly, under the same 
transformation J^(<S) = c [ e t ( -j +At - ) z/ ('S>)? and for A = I, it can be written as u'(S)(t,T) = u(S)/(l + 

t) 2 . 

In any case, the derivative of v(S) w.r.t. matrix A is exactly equal to —\v(S). We also know that 
if the pattern is optimal w.r.t. linear transformation, the derivative of v(S)o~o w.r.t. to matrix A 
shall be null. This implies that v(S)D — lv(S)cro = , which leads to D = ctqI and T = 0. □ 
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Figure 1: Square, Hexagonal and Triangular grids. 



We know that T is symmetric and T = 0. Thus, tr(T) = 0, i.e., Eigen values are invariant by 
rotation. When a grid is optimal, we must have T = and the matrix T must be invariant w.r.t. 
isometric symmetries of the grid. On 2D plane, the grid patterns which satisfy this condition are 
square, hexagonal and triangular grids only. The square grid is symmetric w.r.t. any horizontal 
or vertical axes of the grid and, in particular, with rotation of it/2 represented by J. Therefore, 
the Eigen system must be invariant by rotation of tt/2. This implies that the Eigen values are the 
same and therefore null since tr(T) = 0. Same argument also applies for the hexagonal grid with 
the invariance for 7r/3 rotation and for the triangular pattern with invariance for 27r/3 rotation. 



5 Local Capacity of Grid Pattern Schemes 

We consider that the set S is a set of points arranged in a grid pattern and, for every slot, the grid 
pattern is the same modulo a translation. We have covered grid layouts of square, hexagonal and 
triangle which are also shown in Fig. [TJ Grids are constructed from d which defines the minimum 
distance in-between neighboring emitters and can be derived from the hop-distance parameter of 
a typical TDMA-based protocol. The density of grid points, A, depends on d. However, the local 
capacity is independent of the value of d or A as it is invariant for any homothetic transformation 
of the set of emitters. 

Our aim is to compute the size of the reception area, A{(\, /3,a), around each emitter %. By 
consequence of the regular grid pattern, all reception areas are the same modulo a translation (and 
a rotation for the hexagonal pattern), and their surface area size, <r(A, 0, a), is the same. If Cj(/3, a) 
is the closed curve that forms the boundary of Ai(X, /3, a) and z is a point on Cj(/3, a), we have 

a(X, f3, a) = i J det(z - z h dz) , (5) 

where det(a, b) is the determinant of vectors a and b and dz is the vector tangent to Cj(/3, a) at 
point z. det(z — z^dz) is the cross product of vectors (z — Zi) and dz and gives the area of the 
parallelogram formed by these two vectors. Equation (jSJ) remains true if z\ is replaced by any 
interior point of Aj(A, (3, a). The SIR S{(z) of emitter i at point z is given by ([I]). We assume that 
at point z, Si(z) = (3. On point z we can also define the gradient of Si(z), VSi(z). VSi(z) is inward 
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dz — J| V5 *| 2 || At 

z 



Figure 2: Computation of the reception area of emitter i. 
normal to the curve Ci([3,a) and points towards Zj. The vector dz is co-linear with J .^j 2 ?, where 



J is the anti-clockwise rotation of 3it /2 (or clockwise rotation of 7r/2) given by J 

\VSi(z)\ 



\VSi(z)\ 

1 
-1 



Therefore, we can fix dz = J^^At and in © 



det(z - Zi, dz) = (z - Zi) x J r^Mr , At 

= -(* - «).^U , 
where At is assumed to be a small step size. The sequence of points z(k) computed as 

Z(0) = 2 

gives a discretized representation of Ci({3,a). Therefore, ((5} reduces to 

, ( A, ftQ )»-iB^)-^^f||||A', (6) 

assuming that we stop the sequence z(/c) when it loops back on or close to the point z. Figure [2] is 
the figurative representation of the computation of the reception area of an emitter. 

The point, z(0) = z, can be found using Newton's method. First approximate value of z, required 
by Newton's method, can be computed assuming only one interferer nearest to the emitter i. The 
negative sign in (jSJ) is automatically negated by the dot product of vectors (z(k)-Zi) and VSi(z(k)). 

The local capacity is given by 

c(z,P,a) = E(N(z,P,a)) = N(z,(3,a) = Aa(A,/3,a) , 
where a(A,/3,a) is computed using above described method. 

6 Local Capacity of Slotted ALOHA Scheme 

In slotted ALOHA scheme, the set of simultaneous emitters, in each slot, can be given by a uniform 
Poisson distribution of mean A emitters per unit square area [UEMT3]. Therefore, using the results 
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u x 


Uy 


U X x 


Uxy 

Uyx 


Uyy 


\H(U)\ 


Square 








-.0052 





-.0052 


.00002704 


Hexagonal 








-.0102 





-.0102 


.00010404 


Triangular 








-.0041 





-.0041 


0.0000168 



Table 1: Numerical differentiation (first and second order partial derivatives) of local capacity of grid pattern 
schemes. f3 = 10.0 and a = 4.0. 

from [4J, we can derive the analytical expression for the local capacity with slotted ALOHA scheme. 
The average size of the reception area around an arbitrary emitter is 

a(A,0,a) = i5#^-S. (7) 

A —71 

a 

Therefore, the analytical expressions (j3J) and ([7j) lead to 

c(z, /3, a) — Xa(X, (3, a) — a(l, /3, a) . (8) 

7 Evaluation 

In order to approach an infinite map, we perform numerical simulations in a very large network 
spread over 2D square map of 10000 x 10000 square meters. The emitters are spread over this 
network area in square, hexagonal or triangular pattern. For all grid patterns, d is set equal to 25 
meters although it will have no effect on the validity of our conclusions as local capacity, c(z, (3, a), 
is independent of the value of d or A. 

7.1 Optimality of Grid Pattern Schemes 

The area, A, over which the function U, in (T4j), is evaluated is in the center of the network area 
and is equal to 2500 x 2500 square meters with Zi located at its center: Zi = (xi,yi) = (0,0). We 
set P = 10.0, a = 4.0 and Ax = Ay = 0.1, such that the perturbation in the location of Zi has an 
infinitesimal effect on the reception areas of the emitters on the borders of A. Note that the function 
U is numerically evaluated using the analytical method of $5j The first and second order partial 
derivatives of U are computed according to §4.31 and results, shown in Table [U show that square, 
hexagonal and triangular grid patterns are locally maximum. Similar results can also be obtained 
for all values of /3 and a with the same conclusions. 

7.2 Local Capacity of grid patterns and slotted ALOHA 

In this case, we compute the size of the reception area of transmitter i, located in the center of the 
network ) = (0, 0). The network area is large enough so that the reception area of 
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Figure 3: Local capacity, c(z, (3, a), when (3 is varying and a is fixed at 4.0. 
1 
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Figure 4: Local capacity, c(z, (3, a), when (3 is fixed at 10.0 and a is varying. 

transmitter i is close to its reception area in an infinite map. A depends on the type of grid and 
it is computed from the total number of transmitters spreading over the network area. In case of 
slotted ALOHA, c(z, (3, a) is computed from analytic expressions (J7|) and (IE]). 

Figures [3] and @] show c(z, (3, a) of grid patterns and slotted ALOHA. As a approaches infinity, 
reception area around each transmitter turns to be a Voronoi cell with an average size equal to 
1/A. Therefore, as a approaches infinity, c(z, (3, a) approaches 1. For slotted ALOHA scheme, (JTj) 
and (E]) arrive at the same result. For grid patterns, we computed c(z, (3, a) with a increasing up 
to 100 and from the results, we can observe that asymptotically, as a approaches infinity, c(z, (3, a) 
approaching 1 is true for all protocols. We can also see that the maximum capacity in wireless 
networks can be obtained with triangular grid pattern and its local capacity is at most double the 
capacity of slotted ALOHA. 





Triangular grid - 
Square grid *- 
Honeycomb grid - 
Slotted ALOHA - 
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8 Conclusions 



Our analysis shows that transmission scheme, based on triangular grid pattern is locally optimal. 
Moreover, compared to slotted ALOHA, which does not use any significant protocol overhead, 
triangular grid pattern can only increase the local capacity by a factor of 2. The conclusion of this 
work is that improvements above ALOHA are limited in performance and may have significantly 
higher protocol overheads. Note that, considering the associated protocol overheads, our results 
may encourage network protocol designers to concentrate on designing robust protocols based on 
simpler medium access schemes. In the continuation of this work, we extend this analysis to include 
more practical schemes like node coloring and carrier sense based schemes and see how do they 
compare with the schemes discussed in this article. 
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